Vortex-lattice melting in magnesium diboride in terms of the elastic theory 
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In the framework of elastic theory, we study the vortex-lattice melting transitions in magne- 
sium diboride for magnetic fields both parallel and perpendicular to the anisotropy axis. Using the 
parameters from experiments, the vortex-lattice melting lines in the H-T diagram are located sys- 
tematically by various groups of Lindemann numbers and the anisotropic parameters. It is observed 
that the theoretical result for the vortex melting with parallel and perpendicular fields agrees well 
with the experimental data. 
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I. INTRODUCTION 

Since the discovery of superconductivity with 
39K in magnesium diboride (MgB2) was announced 
it has caused a large number of experiment al H, ^a, 
H, S 0, Q and theoretical investigations [1, llCll |. This 
introduced a new, simple (three atoms per unit cell) bi- 
nary intermetallic superconductor with a record high (by 
nearly a factor of 2) superconducting transition tempera- 
ture for a nonoxide and non-C60-based compound. It dis- 
plays a variety of unusual properties. The high transition 
temperature seems to be either above or at the limit sug- 
gested theoretically several decades ago for BCS, phonon- 
mediated superconductivity jll| . It is characterized by 
a double band superconductor with two superconducting 
gaps of different size, the larger one originating from a 
quasi-2D cr-band and the smaller one from a 3D 7r-band. 
The electronic cr-states are confined to the boron planes 
and couple strongly to the in-plane vibration of the boron 
atoms. The unfixed anisotropic parameter 7 vary widely 
, ranging from 1.1 to 6 dependin g o n the measurement 
technique and on sample types [1, [l^, [H, [13] ■ Although, 
some experimental data can shed light on the mechanism 
of superconductivity, it was demonstrated 15| that the 
type-II superconductivity may not be well described by 
the standard anisotropic Ginzburg-Landau theory, just 
due to the two-band structure. 

The vortex-lattice solid (glass with random pinning) 
state with zero linear resistivity is crucial for the appli- 
cation of high- Tc superconductors, thus the melting of 
vortex-lattice in bulk type-II superconductors is of great 
significance [H, HI [11, Hi, HO] . The traditional Linde- 
mann theory suggests that the lattice melts when the 
root mean square thermal displacements of the compo- 
nents of a lattice reach a certain fraction of the equi- 
librium lattice spacing. This criterion was first adopted 
to study the vortex-lattice melting transition in type- 
II superconductors with a magnetic field parallel to the 
anisotropic axis 17[ , then it was used to draw the melt- 
ing lines when mag netic field is perpendicular to the 
anisotropic axis [2ll . [2^ . [23l [23 | . Since the upper criti- 



cal field is also high in MgB2, the thermal fluctuation 
may drive the vortex-lattice to a vortex liquid in a field 
far below the upper critical field [2l| through vortex melt- 
ing. 

In this paper, using the parameters measured in recent 
two typical experiments 0,13, we study the vortex-lattice 
melting transitions in the framework of the elastic theory 
phenomenologically. The melting lines for magnetic fields 
both parallel and perpendicular to the anisotropic axis 
are systematically located with different groups of Linde- 
mann numbers and anisotropic parameter. The compar- 
ison with experiments are also made. The present paper 
is organized as follows. In Section 2, we introduce the 
theoretical method. Section 3 presents the main results. 
We give a short summary in the final section. 



II. ELASTIC THEORY 

For ideal triangular vortex-line lattice, the free en- 
ergy in elastic theory for fields parallel and perpendicular 
to the X axis can be generally written with quadratic 
terms of the deviation vector u = (ux,Uy) describ- 
ing the fluctuations of vortices from their equilibrium 

positions [H [13, M, [13, [H, m [11 
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The matrix C for parallel fields is different from that for 
perpendicular fields. We denote and C'**' to be the 
elastic matrix for the fields parallel and perpendicular to 
c-axis, respectively, which are given as follows 
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where k"]^ = + ky, cq%,Cl = cn - cge, and C44 are 
the wave-vector-dependent shear, bulk, and tilt elastic 
moduli. The detailed expressions for these moduli can be 
found in several previous papers [l^, [13, [12, [13, [1E| • The 
thermal fluctuations of the vortices are given by inverting 
the kernel C(k) 

(««)-||^ /rfk[CL]-i(k),a = x,y,/3 = c,a5. (4) 

The Lindemann criterion presumes that the lattice 
melts, when the root mean square thermal displacements 
of the components of a lattice reach some fraction of the 
equilibrium lattice spacing. For parallel fields, we con- 
sider the mean-square displacement of a vortex lattice 
from the equilibrium dPiT) = u^ + Uy, then use the usual 
isotropic Lindemann criterion , 

(d^) =c^a\ (5) 

For perpendicular fields, an anisotropic Lindemann cri- 
terion should be employed 

/ 2\ 2 2 / 2\ 2 2 /n\ 

{Uj = C^a^, [Uy) = CyUy. (6) 

with Cx and Cy are two Lindemann numbers for two trans- 
verse directions. Combining the Lindemann criterion and 
the elastic theory, the melting equations are derived with 
the Lindemann coefficients c's. 

If fields are applied perpendicular to c-axis, the effect 
of layer pinning reduces fluctuations in both directions 
and induces an additional momentum-independent term 
to the elastic matrix in Eq. (3) such that 



III. RESULTS AND DISCUSSIONS 

A. Vortex melting lines for fields parallel to the 
c-axis 

First, we calculate the melting line of MgB2 thin film 
for fields parallel to the c-axis. The experiments were 
performed by Wen et al. [3| We employ the experimen- 
tal parameters as follows: the zero-field superconduct- 
ing transition temperature Tc = 39K, The Ginzburg- 
Landau parameters k = 26, the upper critical fields 
along c axis H^2{T = 0) = 15T. Since the anisotropy 
coefficient 7 = H^2/Hc2j i^^t well established for 
MgB2 [§, [H, [H, [i3| • We vary the anisotropic parameter 
7 in the range of 1.5 to 4.5 in our calculation. When fields 
parallel to the c-axis, we set Cx = Cy for no anisotropy 
exists here. We draw the melting lines by elastic theory 
with the Lindemann number c^, = c„ = 0.1 and with the 




T(K) 

FIG. 1: The comparison of vortex-lattice melting lines for 
magnetic fields parallel to c-axis with experimental results in 
Ref. [3]. The anisotropic parameter 7 varies from 1.5 to 4.5 
with the step 0.5( from top to bottom ). The inset shows 
the melting line with the anisotropic parameter 7 = 3 for 
Lindemann number Cx = 0.1. 
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FIG. 2: The vortex-lattice melting lines of MgB'z single crys- 
tal for magnetic fields parallel to c-axis. The Lindemann num- 
ber Cx = 0.1 and the anisotropic parameter 7 varies from 1.5 
to 4.5 with the step 0.5( from top to bottom ). Sohd circles 
are experimental results from Ref. [3|. The inset shows the 
comparison of melting line with 7 = 3 and Cx = 0.1 with the 
experimental results Q|. 



variation of 7, as shown in Fig.l. Interesting, in the inset 
of Fig.l we find that the melting line with 7 = 8 agrees 
well with the irreversibility lines measured by SQUID in 
ref. [3|. The irreversibility line in superconductor is usu- 
ally regarded as the melting line. [2l| . The value of 7 = 3 
in this sample is similar to that on transport measure- 
ments of the upper critical field anisotropy performed on 
single crystals [H, [13] . 

Then, we evaluate the melting line of MgB2 single 
crystal for fields parallel to the c-axis. This experiment 



'P \^ ciikxky ciikl + cl^kl + cl^k'j 

(7) 

where @ is proportional to the critical depinning current 
and depends on the layer separation s [28| . 
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was done by Eltsev et al. [8|LThe experimental param- 
eters measured from Ref. [8[ are collected as follows: 
Tc ^ 38. 5K, K ^ 26, H^^iT = 0) = 7.5T. We also change 
the anisotropic parameter 7 from 1.5 to 4.5 as above. It 
is found that the variation of 7 also has influence on the 
melting line. Fig. 2 shows the theoretical melting results 
with Cx — Cy = 0.1, where the experimental data are also 
exhibited. As shown in the inset of Fig. 2, we observe 
that the melting line with 7 = 3 agrees with experimental 
data 3. Here the resulted anisotropic parameter 7 = 3 
also matches up to the experiment results on transport 
measurements of MgB2 single crystals [H, 



B. Vortex melting lines for fields perpendicular to 
the c-cixis 
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FIG. 3: The comparison of vortex-lattice melting lines for 
magnetic fields perpendicular to c-axis with experimental re- 
sults in Ref. Q (without intrinsic layer pinning). 



We next turn to the vortex melting transition for fields 
perpendicular to c-axis. The experimental parameters for 
a MgB2 single crystal are taken from Ref. [1] as the fol- 
lowing: Tc = 38.5K, 7 = 3, K = 26, H^^{T = 0) = 22r. 
Following the procedure outlined in the preceding sec- 
tion, we calculate the thermal fluctuations in the two 
transverse directions. Setting Cx = Cy = 0.14, we ob- 
tain two curves similar to those in Ref. l22l| . which are 
shown in Fig. 3. As pointed out in Ref. [2J|, to inter- 
pret these two curves as two " melting lines " and thus 
reach to a conclusion of an intermediate smectic phase 
is unphysical, since the elastic theory can give at best 
one single melting line. To impose the same Lindemann 
number along the two directions is of no physical basis. 
In order to achieve a single melting line, we tune the Lin- 
demann number Cj, at 0.1. A good collapse of the melting 
lines in two directions can be achieved if setting the ra- 
tio Cx/cy « 1.4, as shown in Fig. 3. It is interesting to 
note that this ratio is very close to 1.37 observed in Ref. 
[13] using parameters in cuprate superconductors. More 
importantly, we reproduce reasonably the experimental 
melting line in Ref. Q as shown in fig. 3. 

The vortex melting for fields perpendicular to the c- 
axis is also influenced by the layer pinning. In order 
to study this intrinsic layer pinning effect, the matrix 
(7) is used to calculate the thermal fluctuations along 
two transverse directions. Here, the layer separation of 
magnesium diboride is s = 3.524AQ. The melting lines 
for 7 = 3 with Lindemann number Cx = Cy = 0.12 are 
collected in Fig. 4. We flnd that the melting line with 
Cx = 0.12 matches the experimental data. By tuning 
the Lindemann number Cy, a good collapse can also be 
obtained by setting the same ration Cx/cy ~ 1.33. 



IV. SUMMARY 

In the framework of the elastic theory, we have studied 
the melting transition in magnesium diboride for mag- 
netic fields parallel and perpendicular to the anisotropy 
axis. Using the parameters from experiments, the melt- 
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FIG. 4: The comparison of vortex-lattice melting lines for 
magnetic fieldsperpendicular to c-axis with experimental re- 
sults in Ref. [8| (with intrinsic layer pinning). 



ing lines with various Lindemann numbers and the 
anisotropic parameter are composed. Although the 
anisotropic parameter varies with the different measure- 
ment technique and sample types, the melting lines de- 
rived here agree well with those in thin films Q and sin- 
gle crystal[8| by using Lindemann number Cx =0.1 and 
anisotropic parameter 7 = 3 for magnetic field parallel 
to the c— axis. For fields perpendicular to the c-axis, it is 
observed that thermal fluctuations normalized by vortex 
separations in the two transverse directions are propor- 
tional to each other, similar to those observed in cuprate 
superconductors. The ratio Cx/cy to achieve a single 
melting line is very close to that observed in Ref. [24| . 
More interestingly, by using Cx — 0.12 and Cy — 0.09, 
we are able to draw the melting line which fits the ex- 
perimental melting line in Q quite well. Although the 
standard anisotropic Ginzburg-Landau theory may not 
be applicable to magnesium diboride, the elastic theory 
of vortex matter can provide a good description of the 
vortex melting in this two-band superconductors. 
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